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We report new chemoconvective pattern formation phenomena observed in a two-layer system of 
miscible fluids filling a vertical Hele-Shaw cell. We show both experimentally and theoretically that 
the concentration-dependent diffusion coupled with the frontal acid-base neutralization can give 
rise to formation of the local unstable zone low in density resulting in a perfectly regular cell-type 
convective pattern. The described effect gives an example of yet another powerful mechanism which 
allows the reaction-diffusion processes to govern the flow of reacting fluids under gravity condition. 
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In last decades the interaction between reaction- 
diffusion phenomena and pure hydrodynamic instabilities 
has attracted increasing interests both from the funda¬ 
mental point of view of nonlinear science and from the 
chemical engineering m- The interest arises from the 
fact that the chemically-induced changes of fluid prop¬ 
erties such as density, viscosity, thermal conductivity or 
surface tension may result in the instabilities, which ex¬ 
hibit a large variety of convective patterns. 

Scenario for instability development essentially differs 
for immiscible and miscible systems of liquids. The sim¬ 
ple, irreversible chemical scheme such as a neutralization 
reaction A-fB —>8 occurring in binary liquid-liquid im¬ 
miscible systems was studied in [iHa]. The pattern for¬ 
mation in the form of irregular plumes and fingers was 
shown to originate from the coupling between different 
gravity-dependent hydrodynamic instabilities occurring 
when an organic solvent containing an acid A is in con¬ 
tact with an aqueous solution of an inorganic base B [4] . 
This irregularity looked natural since the configuration 
of more dense acid on top of less dense base in the pres¬ 
ence of gravity is unstable via the Rayleigh-Taylor (RT) 
mechanism [SJE]. A much more impressive, regular pat¬ 
tern of cellular-like fingers keeping contact with the inter¬ 
face was reported for an organic base [7]. The regularity 
was shown to originate from the perfect balance between 
the RT instability on the one hand, and the Rayleigh- 
Benard [S] and Marangoni mechanisms on the other. 
Thus, a liquid-liquid interface was recognized to be im¬ 
portant for performing fine-tuning of salt fingers. 

A completely different situation was observed in the 
miscible case. The main engine breaking the equilibrium 
here was found to be the difference between the diffusion 
rates of all three substances resulting in double diffusive 
(DD) instability or diffusive-layer convection (DEC) as 
well as RT instability m- All the works devoted to this 
subject at a given moment, usually noted the formation 
of irregular patterns of fingers. For example, recently 
Almarcha et al. niHn] have shown that the various pos¬ 
sible convective regimes can be triggered by acid-base 
reactions when a less dense acid solution lies on top of a 
denser alkaline one in the gravity field. The possible dy¬ 


namics are a composition of only two asymptotic cases: 
irregular plumes induced by a local RT instability above 
the reaction zone and irregular fingering in the lower so¬ 
lution induced by differential diffusive effects. 

Absolutely in all works in this field, cited or not, the 
diffusion coefficients of species have been assumed to be 
constant. Generally, a concentration-dependence exists 
in most systems, but often, e.g. in dilute solutions, the 
dependence is weak and the diffusion coefficient can be 
assumed constant [14]. This is especially true for fluid 
mechanics m- Some rare examples of the influence 
of concentration-dependent diffusion include the colloid 
ultrafiltration m and membrane transport m where 
the basic fluid flow is just slightly modified. Reaction- 
diffusion problems include the plasma wave dynamics [ 18 j 
and the Turing instability under centrifugal forces [T9|. 

In this Letter, we also focus on the study of chemo- 
hydrodynamic processes which accompany a frontal neu¬ 
tralization reaction taking place between two misci¬ 
ble liquids. We report a new type of instability, the 
concentration-dependent diffusion (hereinafter ODD) in¬ 
stability from the family of the double-diffusion phe¬ 
nomena [20]. It arises when the diffusion coefficients of 
species depend on their concentrations. We demonstrate 
both experimentally and theoretically that chemically- 
induced changes of reagent concentrations coupled with 
concentration-dependent diffusion can produce spatially 
localized zone with unstable density stratification that 
under gravity gives rise to the development of perfectly 
periodic convective structure even in a miscible system. 

Experimental results. The experiments were performed 
in a vertically oriented Hele-Shaw cell made of two glass 
plates (width 2.5 cm x height 9.0 cm) separated by a thin 
gap of 0.12 cm. The cell was filled with aqueous solutions 
of reagents whose concentration always provided a steady 
stratified density distribution. We examined a few acid- 
base pairs formed by HCl or HNO 3 from one side and 
NaOH or KOH from another. During the filling of the cell 
with the upper solution, the lower layer was separated by 
a thin plastic slide inserted in two narrow (0.3 mm) slots 
made in the walls. Fizeau interferometry was used to 
visualize a refractive index distribution. The latter was 
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FIG. 1: (a-c) Chemo-convective structures arising due to the ODD instability observed 1100 s after as the aqueous solutions 

of HNOs (top) and NaOH (bottom) were brought into contact in a vertical Hele-Shaw cell: (a) Velocity field revealed by the 
tracks of light-scattering particles; (b) Interferogram showing a refractive index distribution; (c) pH distribution obtained in 
the presence of color indicator. Initial concentrations of species are both equal to 1 mol/1. The initial contact line is indicated 
by the horizontal band. (d,e) Evolution of the npper (indicated in (a) as Zu) and lower (Zj) boundaries (d) and wavelength (e) 
of the CDD pattern in time obtained experimentally (points) and numerically within the theoretical model (lines). 


caused by inhomogeneities induced by the concentration 
distribution of species and reaction exothermicity. The 
maximum temperature increase measured in the vicinity 
of the reaction front was found to be near 1 K. In this case 
the refractive index deviation due to temperature was at 
least one order of magnitude smaller than that caused 
by concentration. Thus, the interferograms obtained in 
the experiments reflected mainly the concentration dis¬ 
tribution. Silver-coated hollow glass spheres were added 
to the liquids to observe the convective patterns which 
form during the reaction. We visualized also the pH dis¬ 
tribution by adding a small amount of universal acid-base 
indicator to the solutions. The comparison of the results 
obtained with and without the indicator has shown that 
the presence of the indicator did not influence the insta¬ 
bility scenario and pattern formation process, as it was 
demonstrated in some recent studies [2TJ [22] . 

Right after the prepared solutions were brought into 
contact, the transition zone started to form between them 
where the reagents were transported towards the reaction 
front only via the diffusion mechanism. Then the occur¬ 
rence of a depleted layer just above the diffusion zone 
has given rise to the formation of plumes which result 
in the development of weak buoyancy-driven convection 
in the whole upper layer. A few minutes after the be¬ 
ginning of the experiment, the fluid flow in the form of 
the periodic array of convective cells (Fig. has been 
formed within the diffusion zone just above the reaction 
front. The cells were arranged between two areas of im¬ 


mobile fluid which definitely indicated the formation of 
a local “pocket” with the unstable density stratification 
(Fig. [^). One can note that the cellular structure did 
not interact with the convection in the upper layer. The 
observations of pH distribution (Fig. &) have shown that 
above and below the cells band the medium has almost 
homogeneous pH (acidic or alkaline). The pH within the 
cellular pattern is more neutral indicating the accumula¬ 
tion of the reaction product in this zone. Clearly defined 
zone of the intermediate acidity suggests that we deal 
here with some kind of a cooperative phenomenon. 

We found that the structure can exist for several hours 
with the band slowly widening with time (Fig. EJi) that 
results in the wavelength growth (Fig. [^). It is impor¬ 
tant to note that this chemo-convective regime was found 
also in all pairs of reactants used in experiments, but only 
at certain ratio of initial concentrations (Table 1). 

Theoretical model. To describe the observed phe¬ 
nomenon of the CDD instability consider two miscible 
fluids filling a close parallelepiped sufficiently squeezed 
along one horizontal direction to use a Hele-Shaw ap¬ 
proximation [6]. The upper and lower layer are aqueous 


TABLE I: Initial concentrations ratio of different acid-base 
pairs for which the cellular chemostructure has been observed. 


Acid:Base 

liN'o'3:NaOH 

HN03:K0H 

HChNaOH 


1 : 1 

1 ; 1.3 

1 : 0.7 





























3 


solutions of acid A and base B respectively. Right af¬ 
ter the process starts, the acid and base diffuse into each 
other and are neutralized with the formation of salt S 
with the rate K. The system geometry is given by two- 
dimensional domain with a;-axis directed horizontally and 
z-axis anti-directed to gravity. We choose the follow¬ 
ing characteristic scales: length - the gap-width h, time 
- h^jDao^ velocity - Dao/h, pressure - povDao/h^ and 
concentration - Aq. Here Poi define constant 

acid diffusivity, solvent density, kinematic viscosity and 
initial acid concentration respectively. The mathematical 
model we develop consists in the set of reaction-diffusion- 
convection equations coupled to Navier-Stokes equation, 
written in the dimensionless form: 


$ = ( 1 ) 


Sc 


dfA 


6 


J(5', $) = - 12$ - Rad^A - 


-Rbd,B - Rsd,S, 
dtA + J($, A) = VDa{A)VA - aAB, 
dtB + J($, B) = VDbiB)VB - aAB, 


( 2 ) 

( 3 ) 

( 4 ) 


dtS + J{'S,S)=VDs{S)\7S + aAB, (5) 

where J stands for the Jacobian determinant J{F,P) = 
dzFdxP — dxFdzP- Here we use a two-field formulation 
for motion equation, and introduce the stream function 
$ and vorticity $ defined by Q. Eq. ([^ differs from a 
standard Navier-Stokes equation by the additional term 
12 $ which is responsible for the average friction force 
due to the presence of the side-walls. Diffusion terms in 
Eq. (3][5 ) have been written in the most general form mi. 

The problem parameters are the Schmidt number Sc = 
u/Dao, the Damkohler number a = KA^h^/Dao and 
the set of solutal Rayleigh numbers for species Ri = 
glSiA^h^/vDao, i = {a, 6 , s}. Their values for the pair 
HNO 3 / NaOH (see Table 1) have been estimated as fol¬ 
lows: Sc = 10^, a = 10^, Ra = 1.5 X 10^, Rb = 1.8 x 10^, 
Rs = 2.4 X 10^ (see Supplementary information B). 

The boundary conditions for Eqs. (T][5) are 


$ = 0 , di'^ = 0 , ^^A = 0 , d^B = 0 , diS = 0 , ( 6 ) 

where i = {x, z} for side-walls and horizontal boundaries 
respectively. The initial conditions at t = 0 are 

z<0: $ = 0, 9;,$ = 0, A = 0, B = 1; 

z>0: $ = 0, = 0, H = 1, B = 0. (7) 

We found that a concentration-dependence of diffusion 
plays an important role in the pattern formation. Thus, 
in Eq. (3][5) the diffusion coefficients have been assumed 
to be not constant, but depend on their own concentra¬ 
tion: Da{A), Db{B) and Ds{S). In order to evaluate the 
diffusion formulas for the pair HNOa/NaOH, we have 
brought together all the known to us experimental data 




p, density k, wave number 


FIG. 2: (a) Instantaneous base state profiles of the total 

density for the case of constant diffusion (circles) and 
concentration-dependent diffusion (squares) at t = 5; (b) Neu¬ 
tral curves for DLC and CDD instabilities which arise in two 
zones low in density shown in Fig. 


(see Supplementary information A) and have obtained 

Da{A) « 0.158A-b 0.881, 

DbiB) « -0.087B -b 0.594, ( 8 ) 

DsiS) « -0.284S'-b 0.478. 

It follows from Q that the salt is most immobile, com¬ 
pared with acid and base. In addition, the diffusivity of 
salt decreases with the growth of the salt concentration. 
All these factors together produce an interesting effect. 
In order to describe it in term of the buoyancy it is con¬ 
venient to introduce the total dimensionless density: 

p{x, z) = RaA{x, z) + RbB{x, z) + RgS{x, z). (9) 

The base state profiles for the density are shown in 
Fig.§L for two different diffusion laws: (i) the constant 
diffusion with standard table values for coefficients; (ii) 
the concentration-dependent diffusion defined by (|^ . We 
see that the curve has only one minimum above the reac¬ 
tion front in the case (i) versus two minima (above and 
below the reaction front) in the case (ii). The lower mini¬ 
mum enclosed within the regions with a stable stratifica¬ 
tion occurs because of the progressively slower diffusion 
of the salt resulting in its accumulation in or near the 
reaction front. Both minima enable for the potential de¬ 
velopment of the instability in the presence of gravity. 
The upper minimum is the typical for the DLC instabil¬ 
ity (Fig. squares). The lower local minimum on the 
same curve is much more interesting: since it has ap¬ 
peared exclusively due to the concentration-dependence 
of diffusion, we have named it as the CDD instability. 

A nonsteady spectral amplitude problem has been 
solved by the method suggested originally in [53| and 
developed for chemo-convection in [am. Fig.[|) shows 
the neutral curves for the DLC and CDD instabilities. 
At time t « 0.15 the CDD disturbance with a wave num¬ 
ber k « 4.6 is the first to lose stability. Then more and 
more waves are involved into the instability area. The 
DLC instability arises at t « 0.25, and its critical wave 
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FIG. 3: Stream function (top) and total density (bottom) 
obtai ned b y numerically solving full non-linear set of equa¬ 
tions ( 1|8 1 for time t = 3. The line 2 = 0 corresponds to the 
initial contact line between layers. Non-linear development of 
the ODD instabilitiy below the contact line is clearly seen. 


decreases with growth of its concentration (the CDD ef¬ 
fect), it can progressively make a local minimum in the 
density profile (figuratively, “density pocket”). Finally, 
under gravity condition one can observe the development 
of the localized cellular convection within the bulk of the 
motionless liquid (see the lower zone Fig. and Fig.||). 

After the localized CDD structure occasionally was 
found in the pair HNOa/NaOH, we have tested a number 
of other systems and found a similar patterning there. In 
our opinion, it may indicate that the discovered effect is 
of a general nature and should be taken into account in 
reaction-diffusion-convection problems as another tool to 
organize the movement of the reacting fluids. The CDD 
effect should take its place among other instabilities (DD, 
DDD, DLC) of the family of the double-diffusive phenom¬ 
ena introduced in physics over 60 years ago |20j . 

We wish to thank A. De Wit for stimulating discus¬ 
sions. The work was supported by the Perm Education 
Ministry, RFBR (13-01-00508a, 14-01-9602 IrmraLa) and 
Basic research program of UBRAS (#15-10-1-16). 


number at the very beginning is fc ^ 0.75. Eventually the 
maximum growth rate of disturbances in the instability 
balloon is shifted towards longer wavelengths. 

To see non-linear development of the disturbances, the 
problem has been solved numerically by a finite- 

difference method described [8] in detail. Stream lines 
and total density of the pattern at time t = 3 are pre¬ 
sented in Fig. We found that in both zones low in 
density the convection develops independently (Fig. 
top). The most interesting situation is in the lower area 
where the cellular chemoconvection with a perfectly peri¬ 
odic structure induced by the CDD instability has been 
observed. As in the experiment, the boundaries of the 
structure slowly move apart with time (Fig. [^). At 
t — 3 the pattern wavelength is about 4.2, which is in 
good agreement with the experimental data (Fig. Si)- 

Discussion and closing remarks. - The system of misci¬ 
ble fluids when a given solution is placed above a denser 
solution with the fastest diffusing species in the upper 
layer has been classified in [10] as a typical case of the 
DLC instability. It occurs due to the formation of a de¬ 
pleted zone low in density which develops above the ini¬ 
tial contact line while an accumulation zone where the 
density is maximal is obtained below the line. As a re¬ 
sult, irregular DLC fingering is observed which develops 
on both sides of the initial contact line. In our case we 
also can identify the DLC plumes above initial reaction 
front (see the upper zone in Fig. and Fig. [^. But 
below the contact line we meet a new kind of instability. 
It may occur only in the reactive case when an emerging 
component starts to accumulate near the reaction front. 
If its molecules quickly leave the reaction zone, then it 
has no significant influence on the instability scenario. 
But if the diffusion coefficient of the reaction product 
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